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Abstract 

We study an interacting electron gas in a quantum wire within the Hartree- 
Fock random phase approximation. Vertex corrections to the electron spin 
polarizability due to the electronic exchange interaction are important giving 
rise to spin-density excitations (SDE) with large oscillator strength shifted 
to lower energies with respect to single-particle states. The energy of in- 
tersubband SDE oscillates with the number of subbands occupied and has a 
minimum when a subband energy is close to the chemical potential. Intrasub- 
band SDE have a linear dispersion at small wave-vectors. The corresponding 
sound velocity is reduced with respect to the Fermi velocity due to exchange 
interaction within the occupied subbands and exchange screening caused by 
virtual transitions to upper subbands. For intersubband and intrasubband 
charge-density excitations (CDE) vertex corrections are of less importance. 
For only a single subband occupied the screening of CDE and SDE in the 



Tomonaga-Luttinger model due to virtual transitions to upper subband is 
studied, where the virtual transitions are treated within the Hartree-Fock 
approximation. The calculations are in good qualitative agreement with ex- 
periments. 
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I. INTRODUCTION 

In recent years progress has been made in preparation and spectroscopic measurements 
of semiconductor nanostructures where electrons are free to move in only one direction. The 
parent system in fabrication of such structures is a doped two-dimensional quantum well with 
large (~20meV) separation between subbands, so that electrons reside in the lower subbands. 
By applying a gate voltage it is possible to further confine the two-dimensional electron gas 
into a one-dimensional wire. The confinement gives rise to one-dimensional subbands in the 
single-particle electron energy spectrum with a subband separation of several meV. This is 
a very clean system to study since the mean free path and inelastic scattering lengths may 
be longer than the length of the wire which is typically of the order of some /im. As a 
result interesting quantum effects may appear. At one-dimensional electron density about 
10^ cm^^ more than one subband can be occupied. 

In Raman spectraaEl collective spin and charge density excitations of GaAs quantum 
wires were observed. The dispersion of low-frequency intrasubband SDE and CDE (plas- 
mons) is clearly seen in angular resolved spectra.™ The wave-vector dependence of the 
spin-wave energy was found to be linear. The corresponding sound velocity is an important 
parameter which depends on the strength of the exchange interaction between electrons. 
Theoretical calculations based on the Hartree random-phase approximationaEl have shown 
good qualitative and quantitative agreement with experiments for the intrasubband and 
intersubband charge-density excitations. 

The role of the exchange interaction is crucial for the collective intrasubband and inter- 
subband SDE. The direct long-range Coulomb interaction leads to depolarization shifts of 
single-particle excitations (SPE) and to the appearance of collective plasma modes. Similarly 
the exchange interaction gives rise to red shifts of spin-flip excitation energies. As a result, 
collective SDE with a large oscillator strength split off the continuum of SPE. This has been 
shown in Hartree-Fock RPA calculations in quantum well semiconductor structures.ll2rE3 In 
this paper we present analysis of electron excitations in a GaAs quantum wire within the 



time dependent Hartree-Fock approximation. The confinement potential was chosen to be 
parabohc. We consider a wire of a finite length much larger than the electron localization 
length for electrons confined in a parabolic well. The properties of the system have been 
calculated for a wide range of electron numbers and occupied subbands. 

We are here primarily interested in the collective spin-density and charge-density exci- 
tations that are two-particle processes. The Hartree-Fock approximation is well suited to 
describe the high energy excitations and all virtual intersubband transitions. The HF-RPA 
equations are treated totally consistently in the numerical approach, i.e. by first solving for 
the HF single-particle states and then with this HF basis solving the equation for the two- 
particle propagator including vertex corrections. This is important in order not to violate 
conservation laws, e.g. to satisfy the generalized Kohn's theorem for dipole intersubband 
charge- density excitations. In the dipole approximation, intersubband charge- density excita- 
tions at q=0 are simply given by the confinement energy since the perturbation only couples 
to the center of mass motion which is decoupled from the relative motion in a parabolic 
confinement .Ii3 Therefore the intersubband charge-density excitations at q=0 do not provide 
much information about the system. However the intersubband spin-density excitations de- 
pend on the relative motion of electrons with different spin orientation and their energies 
are not given simply by the confinement energy. Indeed we found SDE energies to be de- 
pendent on the electron density. These energies oscillate with the occupation of subbands 
and have a minimum when a new subband starts to be filled. For large wave vectors the 
Landau damping becomes significant for the intrasubband charge-density and spin-density 
excitations. 

The HF-RPA, being a more or less reasonable tool for studying intersubband excitations, 
is not a good approach for low energy intrasubband SDE and CDE in quantum wires with 
one or only few subbands occupied. As knownt3 even weak interactions between electrons in 
one-dimensional systems change their ground state from the Fermi to the Luttinger liquid. 
On the other hand, it was shownt^O that the Hartree random phase approximation gave 
the same dispersion for the low-energy plasmons in the spinless system as the Tomonaga- 



Luttinger model. One can not expect that in a similar way HF-RPA gives correct results for 
the SDE spectrum at small wave vectors when electron spins are taken into account. The 
Tomonaga-Luttinger model (TLM)li3'llaH3 is more appropriate for this purpose. However, 
this model is restricted to essentially low-energy physics and cannot be applied to quasi- 
one-dimensional systems where the coupling between the subbands is siginficant. Even if 
only one subband is occupied electron virtual transitions to higher nonoccupied subbands 
should be taken into account. These transitions cause the direct and exchange intersubband 
screening which renormalizes the intrasubband excitation energies. The HF-RPA can be 
used for evaluation of these screening effects. We thus combine the two approaches and 
treat the low-energy intrasubband excitations in the Tomonaga-Luttinger model, while their 
intersubband screening was calculated with the use of HF-RPA. We calculate the direct and 
exchange screening as a function of the confinement energy. It is shown that the screening 
effects are rather weak, stronger for SDE than for CDE and of the order of the correction 
due the backscattering caused by exchange scattering from one side of the Fermi surface to 
the other for sufficiently small subband separation (~5meV). In an experimental situation, 
these effects may be just as important as the backscattering effect and they are not easy to 
separate. The screening may be varied by changing the parabolic confinement, but keeping 
the electronic density constant. 

The paper is organized in the following way. Section || describes the model treating 
the ground state and excited states within the Hartree-Fock random phase approximation. 
Our results for the intersubband excitations are presented in section |I|. The intrasubband 



excitations are discussed in section |rV| including the calculation of the intersubband screening 
factors within the Tomonaga-Luttinger model. In addition the appendices show detailed 
evaluations of some crucial Coulomb matrix elements and how the excitation energies in the 
system can be found from a symmetric eigenvalue problem in the HF-RPA. 



II. MODEL 

We consider a strictly two-dimensional electron system lying in the x-y plane. The mo- 
tion in the z direction is neglected since we assume that the actual thickness of the well 
is smaller than all other relevant length scales. We use the Hartree-Fock approximation 
to reduce the many-particle Hamiltonian to a single-particle Hamiltonian for each electron 
in an effective potential approximating the electron-electron interaction. The correspond- 
ing Hartree-Fock single-particle propagator is shown in Fig. [T|. For the charge- density and 
spin-density correlation functions we use the corresponding time-dependent approximation 
describing the self-consistent linear response of the two-dimensional electron gas to an exter- 
nal perturbation (Hartree-Fock random phase approximation). We assume that the ground 
state is paramagnetic. The Hartree-Fock equation for a quasi-particle is then 
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where fa is the Fermi occupation factor for a state with orbital quantum number a (the spin 
is up or down) with eigenenergy ea and m* is the effective mass. The dielectric constant of 
the surrounding medium is noted by k. We assume that the confining potential is parabolic 

V,{y) = ^m*u;y (2) 

giving a subband spacing of hujQ. Periodic boundary condition in the longitudinal direction 
of the wire gives a Bloch-type single-particle wave function 

^„fc(x, y) = -j=e'^''i)nk{y) , (3) 

where n is the subband index and the longitudinal wave-vector is k=p-2TT / L.j. with p&Z 
and Lx is the length of the wire. For noninteracting electrons in a simple one-dimensional 
parabolic potential the transverse single-particle eigenfunctions are independent of the lon- 
gitudinal wave-vector 

6 



My) = n- I =Hn ( r ) ^^P 



y' 



(4) 



211 

where the nth Herinite polynomial is denoted by Hn- The electron is localized in the 
transverse direction within the order of the confinement length defined as lQ='\h/{m*ujQ)Y^'^. 
The eigenenergies corresponding to the eigenstates (^ are 

E^, = hu:,{n + -) + ,^. (5) 

The effective single-particle Hamiltonian corresponding to ([l|) is diagonalized using the wave 
functions of the noninteracting electrons @) as a functional basis and the self-consistent 
solutions are obtained by iteration. 

For inelastic light scattering not close to resonance the Raman intensities in polarized 
and depolarized scattering geometries are proportional to the imaginary parts of the charge- 
density correlation function and spin-density correlation function respectively.^^ The charge- 
density correlation function is 

7 /-oo 

X+(q,co) = --J^ dte^^\\[p{q,t),p{-q,0)]\), (6) 

and the spin-density correlation function is 

? /"OO 

X-{q,u;) = --J^ dte^-'{\[a{q,t),a{-q,0)]\), (7) 

where p{q, t) is the Fourier transform of the charge- density operator and o"(g, t) is the Fourier 
transform of the spin-density operator along the spin quantization axis and (| . . . |) denotes 
the thermodynamic average. In fact there are three different possible spin-density correlation 
functions, however this triplet is degenerate for the nonmagnetic system discussed here. We 
denote the real part of the correlation functions (||) and (0) as R^{q,u) and the imaginary 
part as 5*^(9, cj) 

X^iq,uj) = -R'^iq,uj)-tS^iq,u). (8) 

The real and imaginary parts of retarded correlation functions are related by the Kramers- 
Kronig relations which is used in the numerical evaluation of the real part of the correlation 
function, see Appendix El. 



If treated exactly the spectral function of the charge- density and spin-density correlation 
function should satisfy the longitudinal f-sum rule independent of the Coulomb interaction 

r duJu;S^{q,u;)=7r^N, (9) 

j-oo m* 

where A^ is the number of electrons in the wire. The charge- density operator is 

piq) =J2l dre^'^-^ra{r)Ur) (al^^T + c\^c,^) (10) 

a,b 

and the spin-density operator is 

a{q) = J2l dre''^-''rair)Mr) (al^^T - cliC,^) (H) 

a,b 

expressed in terms of quasi-particle Hartree-Fock states where Ca-i destroys a particle in a 
state with orbital quantum number a and spin up. Therefore the correlation functions (^ 

and (0) may be written as 



xHq, c^) = E i^l^'nb) {c\e-'nd)nt,cM , (12) 

abed 

where the two-particle charge- density and spin-density correlation functions 

^Xcdi^) = naT6T;cTrfT(^) ± ^a'\bV,cldi{^) ± 

^aibUcld^i^) +'n^aibi;cidii^) (13) 

are combinations of the retarded two-particle Green's function defined as 
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Here greek indices contain orbital and spin quantum numbers. The two-particle Green's 
function (|1^) are found in the HF-RPA, which neglects correlation effects. The fermion 
closed loop and all ladder diagrams are included in the polarization operator and the 
Hartree-Fock single-particle Green's function are used for the fermion lines, see Fig. ^.E^E^ 
The Hartree-Fock random phase approximation is equivalent to a time dependent Hartree- 
Fock approximation. Since the time dependent Hartree-Fock approximation is a conserving 
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approximation,EZl the longitudinal f-sum-rule (|^) is in general satisfied and serves as a test 
of the consistency of our numerical procedure. By using the Matsubara techniques (in- 
troducing the complex frequency ihuJn='i2n'Kki,T where n is an integral number, kb is the 
Boltzmann constant and T is the temperature) to find the two-particle Green's function at 
finite temperatures, the matrix equation 

^jsi'^^n) 22 i^j'S;-y5' — Vs-y'-^^S') 'n^al3;Y5' (i^n) , (15) 

7'<5' 

is obtained where n°^(ic(j„) is the unperturbed two-particle Green's function. It is given by 
the single-particle Hartree-Fock Green's functions 

U°s{iuJn) = -^J2 G^{iuJn + iuJn')Gs{iUn>) 

fs - A ^^g^ 



ihujn - (e^ - e^) 
and the Coulomb matrix element is 



e 



Mr)Mr)r,ir')Mr') 



V.r,S0 = - dr dr' ^^^' '^' '^^\ '^"' ' . (17) 



The negative term in the summation (|1^) represents the vertex corrections. In the Hartree 
random phase approximation, this term is neglected. We now introduce the induced "den- 
sity" matrices K^^{q,uj) 

KfM.^) = E(«le^''-1&)na%c.(^) , (18) 

ab 

in order that the correlation functions may be expressed as 

xHq.oj) = Y.KfM.^){c\e-'nd) . (19) 

cd 

The charge-density induced "density" matrix satisfies the equation 

KtMM = . ^V~^" [2{a\e-^^^\hY + 
Tiuj - [ea - th) 

E (2Kb;ad - Vbc-M) KUq^ uj)] , (20) 

cd 
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and similarly for the induced spin-density matrix 

EHc;<«iif.-.(9,w)]. (21) 

cd 

If vertex corrections are neglected the spin-density excitation spectra are simply given by 
single-particle Hartree-Fock energies. In HF-RPA spin-density excitations are shifted from 
the single-particle Hartree-Fock energies due to the exchange interaction. As a result collec- 
tive spin-density excitations may appear in the spectra. The shift from the single-particle 
Hartree-Fock energies is a measure of the strength of the exchange interaction in the system. 
For long wavelength charge- density excitations the direct Coulomb interaction dominates 
giving rise to plasmons, so that the exchange interaction is not so important here. It will 
however give a different coupling of the single-particle excitations and collective excitation 
(Landau damping). The eigenequations (^) and (|21| ) have to be solved numerically. See 
Appendix ^ for a more detailed description of the numerical approach. 

III. INTERSUBBAND EXCITATIONS 

In the conventional backscattering geometry for Raman experimentsEJEl the wavelength 
of the incident and scattered light are much longer than the width of the wire. The dipole 
approximation is therefore valid for intersubband excitations. Far infrared absorption mea- 
surements in the transverse direction provide the same information as polarized Raman 
scattering, which in the case of parabolic confinement is simply that the generalized Kohn's 
theorem is satisfied giving a single peak in the spectra at Kuq. Intersubband spin-density ex- 
citations may supply interesting additional spectroscopic information. In the noninteracting 
case and within the dipole approximation the spectral function is 
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and the static susceptibility is 



^°^^^) = (^)/(^^o)^ (23) 



for both the spin-density correlation function and charge- density correlation function. The 
charge- density correlation function does not depend on the Coulomb interaction according 
to the generalized Kohn's theorem in the dipole limit. Therefore S{qy,uj) and R{qy) (qy^O) 
are invariant with respect to the electronic density for charge- density excitations. 

In the calculation of the Hartree-Fock ground state, the functional basis has been chosen 
large enough in order that further expansion of it or further iteration of the Hartree-Fock 
equations does not result in visual change to the single-particle energy spectra or the electron 
density of the ground state. To attain sufficient accuracy in the calculation of the charge- 
density excitations and the spin-density excitations the size of the functional basis of the 
excited states has been chosen such that further refinement results in the change of the 
location of the excitation peaks smaller than a typical linewidth in experiments hrj^OAmeV. 
For the calculations we employ the usual GaAs parameters, m*=0.067r?7,o, where mo is the 
electron mass, and the dielectric constant n=12A. The calculations have been performed 
for r=1.0 K. 

We have calculated the intersubband spin-density excitations as a function of the electron 
density for two different confining potentials; /iu;o=11.37meV (/o=100A) and /ia;o=7.90meV 
(/o=120 A). The length of the wire in both calculations is Lx=24:00 A. In all our results, the 
generalized Kohn's theorem is satisfied with a high degree of accuracy for charge-density 
excitations, i.e. there is only a single peak in the spectra at huQ having all the spectral 
intensity satisfying the f-sum rule (j^) as given in (^2]). The generalized Kohn's theorem is also 
satisfied within the Hartree random phase approximation. Therefore a Hartree RPA and a 
Hartree-Fock RPA calculation will give the same result for the long wavelength intersubband 
charge- density excitations for a parabolic confined wire. 

Numerical results for the spin-density excitations as a function of the electronic density 
are shown in Fig. ^ where ^C(Jo=11.37 meV and Fig. | where ^0^0=7.90 meV. The upper 
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panel of Figs. |^ and ^ shows the static spin-density susceptibihty and the center graph the 
excitation energy of the dominant peak in the spectra (SDE). The lowest subfigure shows the 
difference in energy between states at the zone center and the chemical potential for different 
subbands (e„,fc=o-/^)- States below the horizontal hne are occupied. The dominant peaks in 
the spectra have intensities that are typically one order of magnitude larger than the weaker 
peaks in the spectra in the density range shown. The intersubband SDE energy is red-shifted 
with respect to the single-particle Hartree-Fock energies. The excitation energy oscillates 
with the subband filling and has local minima where a new subbands starts to be filled. In 
this case the screening is large so that the self-consistent subband spacing is significantly 
smaller than the confinement energy hujQ. In Fig. ^ the spin-density excitation with the 
largest oscillator strength has been calculated for densities common in present experiments 
on wires. In this calculation only, we use the local density approximation (LDA) for the 
exchange interaction with the parametrizationE^ 



V-fe).-J^:^£!. (24) 

V TT K 

For high density the oscillations of the spin-density excitations are weak. When the 2D 
electronic density is sufficiently low and we are close to a new subband starting to be 
filled, we see that the HF-RPA gives overdamped modes (imaginary excitation energy, see 
Appendix |B[) for the spin-density excitation, e.g. for /ici;o= 11-37 meV (Zo=100 A) there is 
an overdamped region ni£)=(7.1-8.3)xl0^ cm~^ and for huJo=7.90m.eV (/o=120 A) there 
are overdamped regions ?t,i/)=(4. 6-7. 5, 12.1-14.5) x 10^ cm^^ for the densities shown. The 
presence of overdamped modes means that our Hartree Fock ground state is not stable in 
these density regimes. The oscillations in the excitation energy are followed by peaks in the 
static paramagnetic susceptibility and the static paramagnetic susceptibility is divergent in 
the overdamped regions. By comparing the results in Fig. (0) and Fig. (^) we see that 
there are broader unstable regimes for /ico'o=11.37 meV (/o=100A) than for /iti;o=7.90 meV 
(/o=120A) as should be expected since the electronic density is lower in the latter case and 
the validity of the HFA is not so good. Almost the same instability regions are found in the 
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case of the LDA where correlation effects are neglected (|2^). 

IV. INTRASUBBAND EXCITATIONS 

When the wave vector component along the wire is not zero, intrasubband excitations 
shows up in the spectra. We will first consider the case of three subbands occupied. For 
the case of two subbands a comparison with experiments has been briefly communicated in 
Ref. ^ Finally we study how the virtual transitions to upper nonoccupied subbands may 
alter the elementary excitation energies in the Tomonaga-Luttinger model, for only a single 
subband occupied. 

A. Low- and high-energy excitations: Numerical results 

The external parabolic potential was set to 11.37 meV (/o=100 A) resulting in a self- 
consistent Hartree-Fock subband spacing between the lowest subbands of 6.6 meV and be- 
tween the next lowest subbands of 8.3meV at the zone center. Three subbands were occupied 
with densities no=9.4xl0^cm~^, 721=7.3 x 10^ cm~^ and n2=3.3xl0^cm~^. We take the wire 
to have a length oi Lx=1.0 /im. 

For three subbands occupied it is generally expectedoQ that the SDK dispersions are 
linear in q, with sound velocities vq, vi and f 2 for intrasubband excitations in the three lowest 
subbands. The Fermi velocities are given by the subband densities Vi=7Thni/m* (2=0,1,2). 
The CDE dispersion also has three branches. The in-phase mode is approximately^ 



<(g) = \q\,j2{vo + Vi + V2)V{q)/h7r, (25) 

where V{q) is the Fourier transform of the Coulomb potential (e^/fi;)[(a;-a;')^+/Q]^^/^. In 
addition there are two out-of-phase modes with linear dispersion.^ 

We note in general good qualitative agreement of our numerical results with the measured 
Raman spectra.0'H At finite wave-vectors along the wire we found low-frequency SDK and 
CDE corresponding to one- dimensional intrasubband motion of the electron gas. The SDE, 
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as can be seen in Fig. ^ follow a linear dispersion in the long wavelength limit. For three 
subbands occupied we see the three collective spin density modes with group velocities at 
vq, Vi and V2- The Landau damping destroys the higher energy modes at large enough wave- 
vectors. The Landau damping shows up in enhanced intensities of satellite single-particle 
peaks around the SDK, as shown in Fig. |^, since we considered a wire of a finite length 
{Lx = 1.0 fim). For infinite wires one should expect a large number of SPE peaks merging 
into a broad band. We found the sound velocity vq larger than the corresponding Fermi 
velocity ('Do=l-02fo), and the sound velocities vi and V2 smaller than the corresponding Fermi 
velocities, (wi=0.98fi and £'2=0. 64^2). This is caused by the exchange interaction within 
each subband and intersubband coupling giving rise to renormalized sound velocities.lHHlc2l 
The other factor comes from intersubband virtual transitions due to exchange interaction 
leading to exchange screening of SDK. These low-energy features within the HF-RPA must 
be considered as approximate results. 

The q dependence of the intrasubband plasmon energy (CDE+) can be fitted by q[- 
Ing]"^/^, as shown in the inset in Fig. |^. That is what is expected for one- dimensional 
wires.Qc3 For three subbands occupied there must also be two CDEs having linear dispersion 
which may also be seen in Fig. 0. The intensities of these excitations are however weak and 
maybe difficult to observe experimentally. The plasmon decays at higher wave-vectors due 
to the Landau damping, as seen in Fig. ^ 

The calculations performed in Refs. |^-|^ within the Hartree random phase approximation 
have been demonstrated to agree very well with experimental results for the charge-density 
excitations. In order to see the effect of the vertex correction on the intrasubband charge- 
density excitations we have also performed a calculation with the same sample parameters 
as above, but omitting the exchange interaction (not shown). Indeed, we found that the 
vertex corrections have negligible effect on the in-phase intrasubband plasmon energy. It 
agrees within the numerical accuracy (1 %). For the two out-of-phase modes the Hartree 
calculation gives a Fermi velocity 15 % smaller for the highest energy out-of-phase mode and 
almost identical for the lowest out-of-phase mode as compared to the HF-RPA calculation. 
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The HF-RPA generally gives about twice as large relative intensity of the out-of-phase modes 
relative to the in-phase intensity as compared to the Hartree random phase approximation. 
It should be noted that according to our analysis of the wave-vector dependence of 
intrasubband SDE the frequency of the lowest energy SDE decreases at large q<k2, and at 
some value of q close to 2/^2 it goes to zero. That means an intrinsic instability of the system 
with respect to formation of spin density waves (SDW) below some critical temperature. 
This is the well-known Peierls instability of one-dimensional systems. OverhauserEj showed 
that the HF paramagnetic state in 3D systems is always unstable with respect to formation 
of a static spin-density wave having a wave vector q^2kF, and MacDonaldBa found that 
the same was expected in 2D systems. Therefore the Hartree-Fock spin susceptibility of 
the paramagnetic state has a singularity near q=2kp. The instability in the quasi ID wire 
emerges from our mean-field analysis, and in the framework of this approach one should 
take into account SDW long-range order in the ground state. However, exact results making 
use of the Luttinger model predict no SDW long-range order, but predict a slowly decaying 
Wigner crystal of short range order at Akpx^ 

B. Screening in the Tomonaga-Luttinger model 

In the Tomonaga-Luttinger model, the Coulomb interaction is restricted to processes 
within the lowest occupied subband and the kinetic energy is approximated by a linear 
dispersion.lH2l In this way the elementary low-energy excitations may be found exactly by 
using a bosonization technique. We will now study how the excitation energies for the 
Tomonaga-Luttinger model may change when the coupling to the higher subbands is taken 
into account. The virtual transitions to the upper subbands will give rise to direct and 
exchange screening. A second order perturbation argument shows that contributions from 
intersubband transitions will be inversely proportional to the subband spacing. We will now 
give some more quantitative results derived from an effective low-energy Hamiltonian where 
the virtual high-energy intersubband excitations are treated within the HF-RPA. In quantum 
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wires used in Raman measurementsoQ, coupling to higher subbands is always present and 
may not be small. Since this coupling may renormalize the excitation energies appearing in 
the experimental system it is an important physical parameter. The renormalization factors 
may be measured directly in an experiment if one is able to control the subband spacing 
(i.e. the curvature of the parabolic confinement) while keeping the ID electronic density of 
the lowest subband constant. 

For the low-energy excitations the interaction may be approximated by 

V= Y. [H''{k,p)cliP-<imcUp + q/2)cl{k + q/2)cosik-q/2) 

ij,qkp,sa 

- F^^{k,p)cl{k + q/2)c,,{p + q/2)cl{p - q/2)c^,{k - q/2)] , (26) 

where the first term represents the Hartree interaction 

H'^{k,p) = V'^'^^ip, k, k,p) (27) 

and the second the Fock interaction 

F'^{k,p) = V°''°^{k,p,k,p) (28) 

between the lowest subband and the higher subbands. The matrix elements are evaluated 
within the Hartree-Fock basis, V^^'"^^{kp; qr)=f dr f dr' ^*^^{r)'ijjn,rir)'ijj*^p{r')ipm,q{r') / \r — 
r'\. By adding this interaction part to the TLM Hamiltonian the virtual transitions to the 
upper subbands are included. Following the standard bosonization procedurecBEa we now 
introduce the boson operators for charge-density excitations 6+ and spin-density excitations 
b~ in the lowest subband (g > 0), 



Y: {i^ik + q/2)c,^{k - q/2) ± cj^(fc + q/2)c,^{k - q/2)) = b^{q)J^ (29) 

fc>o ' '^ 



Y (dl^ik - g/2)cot(A: + q/2) ± c^ik - q/2)c,^{k + q/2)) = b^{-q)J^ . (30) 



and 



fc<o * ^ 

The spin and charge parts of the Hamiltonian decouple so that the interaction Hamiltonian 

may be written as V=V^+V~ , 
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V^=Y. ^t,ik, q) {X^'\k, q)h{q) + X^^K qP{-q)) , (31) 

ij,kq 

where we have introduced the intersubband charge-density and spin-density operators 

^J(fc, q) = QT(fc - g/2)c,t(fc + g/2) ± 

Qi(A;-g/2)c,i(fc + g/2), (32) 

the charge-density matrix element (n=l,2) 



X:Hk,q) = Jl^l [HHkF,k) - \F:^{k,q) ) , (33) 



and the spin-density matrix elements 



X-HkA) = ^\'^-^\{-\Fi^{Kq)). (34) 

Here F^{k, q) are combinations of the Fock matrix at kp and -kp 

Fi\k, q) = F'^{kp, k)e{q) + F'^{-kp, k)e{~q) (35) 

and 

F'^\k, q) = F'^ikp, k)e{-q) + F'^i-kp, k)e{q) . (36) 

There are also spin-flip terms in the interaction Hamiltonian similar to the spin-spin term 
in ([31| ) , but the spin triplet is degenerate so we can limit the discussion to only one of the 
modes. The TLM may be solved by performing the Bogoliubov transformationcy'Ea 

hHq) + h^K-q) = |^ {^Hq) + ^P^\-q)) (37a) 

h^{q) - h^\-q) = M| (^±t(^) _ ^±(_^)) , (37b) 

where uj{q)=vp\q\ is the unperturbed excitation frequency and E^{q) is the excitation 
energy of the charge-density and spin-density excitations given by the one-band Luttinger 
model. In order to find the elementary intrasubband excitations in the system we define the 
time-ordered boson Green's functions 
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D^{q,u;) 



'i\TP^{q,t)P^\q,0)\) 



(38) 



and 



Diiq,u) 



e'-\\TP^\-q,t)f3^^q,0)\) 



±t( 



(39) 



where the ordinary unperturbed (couphng to the higher subbands is neglected) Green's 
functions is DQ{q,uj) = l/{huj-E^{q)) and the anomalous unperturbed Green's function is 
D^Q(g, C(j)=0. The interaction Hamiltonian (|3T|) may be represented in the same form as 



m 



v^ = ll [ft\SHq) + ft\S^K-q) 



(40) 



where fi^{q) and f2'^{q) are given by the intersubband operators (|5^), the matrix elements 
(^) and (34), and the transformation (|37|) . To lowest order in the intersubband coupling 
the Dyson's equations are 



D^iuj) = D^iu) + 



D^{u)\P^^D^{u) + Pt,Di{u) 



(41) 



and 



Di{uj) = Q + D^{-uo) P,2Di{ij) + PnD^{Lj 



(42) 



where the polarizations P^(g,cj) (q^,/3 = 1,2) are defined as 

Pa^iQ,^) = -j-J_Jte^-'{\Tf^{-q,t)fp{q,0)\) (43) 

and may be expressed in terms of the time-ordered two-particle operator H^^ given by 
similar equations as the retarded two-particle Green's functions 11^ in (p!3|). Intersubband 
excitations are nearly vertical and the intersubband excitation energies are much larger 
than the low-energy intrasubband excitation energy, hence we approximate the polarizations 
(^) by it's limit when {q,uj)-^{0,0). The poles of the Green's function ( PB] ) defines the 
renormalized intrasubband excitation energies 



E^{qf = E^{qf + 2P,,E^{q) + (Pf^ " A'l) • (44) 

The virtual transitions to the higher subbands are included by the polarization operator 



(^) that describes the coupling between the lowest occupied subband and all the higher 
subbands. By carrying out this calculation we find that the shift of the excitation energies 
may be represented by the dimensionless parameters (?t,=1,2) 

IT 

ht = ^ E W^'Kk)W^''^{p)IllUKp), (45) 

^ tjlm,kp 



where Wi ^''{k) and W2 ^\k) are given by 

Wt''{k) = Wi{kF,k) 



and 



}^(^r^[kF^k) + F'^{-kF,k)) , (46) 



W,-'\k) = -- {F^^ikp, k) + F'^{-kp, k)) , (47) 



W^'\k) = -i {F'KkF, k) - F'^{-kF, k)) . (48) 



The intrasubband excitation energies are 



2 



'E^il)\ _ (, f^^i<l)^2i±\ (, r± 



\E-(,)I V ' EH,)' "^^^ -''-)' f^'" 

where hf and k^ represent the contribution caused by screening due to the virtual excitations 
to the nonoccupied upper subbands. Stability of the system with respect to compression 
requires the screening factors not to be too strong. If not, then the perturbation approach 
is not valid since we see from (^9|) that too strong screening factors lead to imaginary 
intrasubband excitation energies. For charge- density excitations the plasmon energy has 
the low-energy behavior g[-lng]^/^ so that the renormalization factor h^ is unimportant. 
The screening factors /i2 are given by differences of the Fock matrix elements with respect 



to kp and -kp (^) and are therefore expected to be small. 
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We have evaluated the intersubband screening factors hi and /i2 for a wire of density 
ni£i=2.0xl0^ cm~^ and for subband separation starting at 2 nieV to 46 nieV, as shown in 
Fig. d^). The screening factor /ig appears in the calculation orders of magnitude smaller 
than hi and may be neglected, as already pointed out above. Even at an external potential 
separation of 2 meV they are smaller than 1%. Since the effect of the screening factor hf 
is unimportant, we find that the plasmon energy is rather insensitive to the coupling to 
higher subbands. Screening has a larger effect on spin-density excitations. The effect on 
the SDE sound velocity is about 5% down to a subband spacing of 4 meV. For smaller 
subband spacings, the exchange screening factor hi appears to decrease, but this is due to 
the fact that the HFA is not valid in this small subband separation regime, because the 
spin-density excitations are overdamped. The screening parameters /i^'^decrease inversely 
proportional to the subband spacing as expected, but the factors /i2 decrease faster than 
inversely proportional to the square of the subband spacing. The screening parameters are 
not very sensitive to the ID electronic density. A similar calculation for niD=3xlO^ cm~^ 
(not shown) shows that the factor hf decreases by less than 10% at huJo=5raeV as compared 
to the results in Fig. (^. Similarly the factor hi decreases by less than 20% at hujQ=5raeV. 

V. CONCLUSIONS 

The importance of the vertex corrections (exchange interaction) that give rise to collective 
spin-density excitations in quantum wires is demonstrated. For charge- density excitations 
vertex corrections are less important. The calculated spin-density excitation spectra and 
charge- density excitation spectra seem to agree well with experiments. The dipole inter- 
subband spin-density excitation energy oscillates with the number of occupied subbands 
accompanied by peaks in the static spin-density susceptibility. For a single subband oc- 
cupied screening by virtual transitions to the upper subband renormalizes the spin-density 
energy by the same order of magnitude as the correction due to backscattering mixing right- 
and left-moving modes. 
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APPENDIX A: SINGULAR CUT-OFF 

In an infinite wire the ground state energy of an interacting electron gas without back- 
ground charges is divergent, which requires a cut-off introduced by background charges 
for the low-wave vector interaction q similar to the 3D situation. Consider the Hamilton 
operator in the basis of the noninteracting states 

^„fe(x, y) = ^=e*^>„(i/) . (Al) 

The Coulomb interaction in the noninteracting basis is 

^mt = ^ E E E yr''{q)iAk - q)clAk' + q)ciAk')c,.{k) (A2) 

imlj kk'q era' 

where Cicr{k) destroys an electron in a state with transverse quantum number i, spin a and 
longitudinal wave-vector k and the Coulomb matrix elements are 

rim;lj 



vr'\q) 



e 



-L_ T — T \ 



which in the limit of long wires can be written as 



dy I dy'KoMy - y')\)My)Uy)My')My') ■ (A4) 
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For small arguments the Bessel function is asymptotically Ko{x)^-\nx. For small wave 
vectors q we may therefore take away the divergent part in the matrix element 



,y-y 



I 



For a finite system we may therefore write 



dy / dy'\n\'—^\ct>,{y)(t>,{y)(t>M)My')- (A5) 



V^^-;'i(g) = 5,^o5,,,5^,iC + V'^^-^'^^q) (A6) 



where C is a large constant. This first term in ( |A6|) gives a part C{N'^-N) to the Hamiltonian 
where A^ is the number operator of particles. We consider states where the total number of 
particles is a good quantum number so that this term is cancelled by a similar term from 
the background charges. The renormalized matrix element for finite wave vector q is given 
by (0) and for q=0 

V'""'^^ {q = 0) = X 

dy J dy' \Yi\^—\(j)i{y)<pj{y)<pm{y')<Pi{y') ■ (A7) 



In the Hartree approximation this is the same as the long wire limit in Refs. [ |3^, |3^ . For 



an easier numerical evaluation by using the integral form of the Bessel function 






^'■»<^'=L*VFTT' <*«' 



the integral (|A4|) may be written as 



du =Iij{u)Irn,l{-u) , (A9) 

where the function Iij{u) is the Fourier transform of the product of a pair of transverse 
wave-functions 

oo 

iuvt 



oo 



VH= / rft;0Kt;)0r(^)e*™*, (AlO) 
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and u and v are dimensionless variables. The function ( |A1CI|) may be written as a linear 
combination of Laguerre functions, so only a one- dimensional integral has to be evaluated 
numerically in (|A9| ). 

APPENDIX B: DIAGONALIZATION OF EXCITED STATES 

The eigenequations (^0[) and (|21| ) are not well suited for a numerical diagonalization since 
they are in general not symmetric. In the case of the time dependent Hartree approximation 
this problem is easily solved due to the local nature of the approximation.^^ The nonlocal 
character of the Fock term causes some problems. However the eigenequations may be 
transformed to a symmetric form as will be described since this may be of interest for 
other calculations based on the Hartree-Fock random phase approximation. This appendix 
therefore covers some of the more technical aspects of the calculations. 

By decoupling (|20D or (l2l| ) into symmetric and antisymmetric combinations the eigen- 



value equation may be rewritten. We define the new variables 

{(^a - eb)ifb - fa 



Uab KQ^^J - — ;, iBiaj 



l/a1'^(9, ^) = {K^biq. u;) - Ktb{q. uj)) ^^-^ (Bib) 

and 

hl,{q) = 2^ifb~fa)iea-eb)[ 

(ale^^^'n^) + {b\e-'''-''\a)], (B2a) 

biiq) = 2./AZZ« [{a\e-^^-\b) - {b\e~^na)) (B2b) 

V Ca — Cfe ^ ^ 

and restrict the summation to e^ < €b obtaining the matrix equations 

- Loy^'^iq, co)=S^- y'^'^iq, u) + b^'^iq) (B3a) 

- ivy'^'^iq, u;) = A^- y^'^{q, uo) + 6^'±(g) , (B3b) 
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where the matrices (since all Coulomb matrix elements are real) S^ and A^ are symmetric 
and antisymmetric combinations of the Coulomb interaction elements defined by 

Kb,cd = \/{fb - fa){fd - fc){ea - eb){ec - ed)[ 

'2(ycb;ad + Vca;bd + Vdb;ac + Vda;bc) 
- {Vbc;ad + Vac;bd + VM;ac + Kd;fec)] (B4) 

and 



A. 



ab.cd 



ifb - fa){fd - fc) r 

^ (ca -e6)(ec-ed) 

'^{Vcbiad — Vca-bd — Vdb;ac + Vda-bc) 
— {Vbc;ad — Vac-fid — Vbd-.ac + Kid;fec)] (B5) 

Similarly for the spin-density excitations the Hartree terms are absent. By now assuming 
that the matrix A is positive definite (this assumption has always been satisfied in our 
calculations), we can use a Cholesky decomposition A=C-C^ where £ is a real matrix and 
diagonalize the system. The excitation energies are given by the symmetric eigenvalue 
problem 

C^ ■S-C-Z = Z-ljI, (B6) 

where the matrix Z contains the eigenvectors of the corresponding eigenvalues in the diagonal 
matrix u;^. If at least one eigenvalue in the matrix ul is negative, it means that we have 
overdamped modes. The Hartree-Fock ground state is then unstable. 

The imaginary part of the charge- density and spin-density correlation functions can now 
be found to be 

S{q,uj) = -lm{x{q,uj)) 

= ±jE^(^Tc.:) {w!{q)TW^^{q)y , (B7) 

X 

where the weights W^{q) and W^{q) are defined as 
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W!{q) = E b'y{q)Ly.2,.^l/u* , (B9) 

yz 

and the indices x, y, z denotes electron-hole pair excitations. From ( [B7| ) it is seen that the 
charge- density and spin-density excitation intensities are positive as they should be. In 
general, the spectral function should have the symmetry S{q,u!)=-S{-q-uj) since from (^ 
and d^) we see that x{q, ^)=X* {-q-^) ■ Here the weights have the symmetry W^{q)=W^{-q) 
and W^{q)=-W^{-q) so that (p7|) has the correct symmetry. 
The real part of the correlation function 

R{q,uj) = -Re{x{q,^)) (BIO) 

may be found from the imaginary part by the Kramers-Kronig relation 

1 t°° 1 
R{q,u) = - du'S{q,u')P , (Bll) 

71 J-oo Uj' — UJ 

where P denotes the principal part. In the static case the susceptibility is therefore 

^„)_ij3(]eM)!±MM!. (B12) 

2 X ^x 

The static susceptibility diverges if there are excitations with very low excitation energy and 
finite weight. 
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FIGURES 
FIG. 1. Feynman diagrams for the single-particle Green's function in the Hartree-Fock Ap- 
proximation. The thick lines represent the Hartree-Fock single-particle propagator, the thin lines 
represent the noninteracting single-particle propagator and the dashed line is the electron-electron 
interaction. 

FIG. 2. Feynman diagrams for the two-particle Green's function in the Hartree-Fock ran- 
dom phase approximation. The shaded boxes represent the four point vertex function, the thick 
lines represent the Hartree-Fock single-particle Green's function, and the dashed line is the elec- 
tron-electron interaction. 

FIG. 3. Intersubband spin-density excitations in the dipole approximation as a function of the 
electronic density. The upper panel shows the static susceptibility, the center graph the excitation 
energy of the dominant peak in the spectra, and the lowest subfigure the difference between the 
single-particle energy at the zone center and the chemical potential. The confinement energy is 
na;o=ll-37 meV and the length of the wire L^=2A00 A, r=1.0K, m*=0.067mo and k=12.4. 

FIG. 4. The same as in Fig. (^), except that the confinement energy is hujQ = 7.90 meV. 

FIG. 5. The energy of the strongest intersubband spin-density excitation in the dipole ap- 
proximation as a function of the electron density for hLOo=7.90 meV (upper panel), and fiuJo=11.37 
meV (lower panel). Only here, the local density approximation has been used for the exchange 
interaction. Other parameters are as in Fig. ^. 

FIG. 6. Intrasubband spin-density excitations. Three subbands were occupied in the wire of 
length La;=1.0/im with densities no=9.4xl0^cm~^, ni=7.3xl0^cm~^ and 77-2=3. 3x 10^ cm^^. The 
external parabolic potential was fic<;o=ll-37meV (/o=100A). The inset shows the dispersion of the 
collective excitations. 
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FIG. 7. Intrasubband charge-density excitations. Three subbands were occupied in the wire 
of length Lx=1.0 /um with densities no=9.4xlO^ cm"^, ni=7.3xl0^ cm~^ and n2=3.3xl0^ cm~^. 
The external parabolic potential was /iaJo=ll-37 meV (/o=100 A). The inset shows the dispersion 
of the collective excitations. 

FIG. 8. Screening parameters {h^ , h2 , h^ and /i2^) as a function of the external parabolic 
confinenient energy Tiujq. The dashed lines show the charge-density screening parameters h^ and 
h2 where h^ is larger than /12". The straight lines show the spin-density screening parameters /i^ 
and /i^ where h^ is larger than h2 . The electronic density of the wire is n=2.0xl0^ cm^-*^. 



29 




5 10 15 

niD(10^cm''') + 



> 

LU 






§ 




2.5 5 

E (meV) 



7.5 



